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ABSTRACT 

Long-term monitoring of the radio emission from supernovae with the Very 
Large Array (VLA) shows that the radio "light curves" evolve in a systematic 
fashion with a distinct peak flux density (and thus, in combination with a 
distance, a peak spectral luminosity) at each frequency and a well-defined time 
from explosion to that peak. Studying these two quantities at 6 cm wavelength, 
peak spectral luminosity (L 6 cm pca k) and time after explosion date (to) to reach 
that peak (t 6 cm pca k — to), we find that they appear related. In particular, 
based on two objects, Type lb supernovae may be approximate radio "standard 
candles" with a 6 cm peak luminosity L 6 cm pca k ~ 19.9 x 10 26 erg s" 1 Hz -1 ; also 
based on two objects, Type Ic supernovae may be approximate radio "standard 
candles" with a 6 cm peak luminosity L e cm pcak « 6.5 x 10 26 erg s" 1 Hz -1 ; 
and, based on twelve objects, Type II supernovae appear to obey a relation 
Le cm peak — 5.5 x 10 23 (t 6 cm pcak - t ) 1A erg s" 1 Hz" 1 , with time measured in 
days. If these relations are supported by further observations, they provide 
a means for determining distances to supernovae, and thus to their parent 
galaxies, from purely radio continuum observations. 

With currently available sensitivity of the VLA, it is possible to employ 
these relations for objects further than the Virgo Cluster out to ~ 100 Mpc. 
With planned improvements to the VLA and the possible construction of more 
sensitive radio telescopes, these techniques could be extended to z ~ 1 for some 
classes of bright radio supernovae. 



Subject headings: supernovae, distances 
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1. Introduction 

A series of papers published over the past 15 years on radio supernovae (RSNe; see 
references to Table 2) have established the radio evolution for 16 objects: 2 Type lb 
supernovae (SN 1983N, SN 1984L), 2 Type Ic supernovae (SN 1990B, SN 19941) and 12 
Type II supernovae (SN 1970G, SN 1978K, SN 1979C, SN 1980K, SN 1981K, SN 1985L, 
SN 1986E, SN 1986 J, SN 1987A, SN 1988Z, SN 1993 J, and SN 1996cb). In this extensive 
study of the radio emission from supernovae (SNe), in most cases using the Very Large 
Array (VLA)fj| radio telescope, two effects have been noticed: 1) Type lb and Type Ic SNe 
have roughly constant 6 cm radio luminosities at peak with only slight differences between 
the two types, and 2) Type II SNe appear to have a higher 6 cm radio luminosity at peak if 
they take longer to reach that peak. These two effects are also consistent with the tenets of 
the SN shock/circumstellar medium interaction model of Chevalier (1982a, b; 1984) for the 
origin of the radio emission if: 1) the radio emission for the rather homogeneous classes of 
Type Ib/c SNe arises in similar circumstellar environments, and 2) the radio emission for 
the more inhomogeneous class of Type II SNe arises in diverse circumstellar environments, 
but the denser the circumstellar medium the longer it takes for the material to become 
optically thin to the radio emission and the brighter the radio emission. If these effects 
can be quantified, they can provide a purely radio-based secondary distance indicator for 
SNe and, by association, for their parent galaxies. The RSNe available for study and some 
relevant properties of their parent galaxies are listed in Table 1. 



4 The VLA is a telescope of the National Radio Astronomy Observatory which is operated 
by Associated Universities, Inc., under a cooperative agreement with the National Science 
Foundation. 
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2. Models 

While the relations to be described are empirical, there is some support from 
theoretical modelling. All known RSNe appear to share common properties of: (a) 
nonthermal synchrotron emission with high brightness temperature; (b) a decrease in 
absorption with time, resulting in a smooth, rapid turn-on first at shorter wavelengths 
and later at longer wavelengths; (c) a power-law decline of the emission flux density with 
time at each wavelength after maximum flux density (absorption r ~ 1) is reached at 
that wavelength; and, (d) a final, asymptotic approach of spectral index a to an optically 
thin, nonthermal, constant negative value (Weiler et al. 1986, 1990). Chevalier (1982a, 
b) has proposed that the relativistic electrons and enhanced magnetic field necessary 
for synchrotron radio emission are generated by the outgoing shock wave from the SN 
explosion interacting with a relatively high-density envelope of ionized circumstellar 
material surrounding the presupernova star. This dense cocoon is presumed to arise from 
mass loss in a stellar wind from a red supergiant SN precursor, or a companion, which 
has been ionized and heated by the initial UV/X-ray flash of the SN explosion. This 
circumstellar material (CSM) is also the source of the initial absorption. A rapid rise in the 
observed radio flux density results from the shock overtaking more and more of the wind 
material, leaving progressively less of it along the line of sight to the observer to absorb the 
slowly decreasing synchrotron emission from the shock region. 

The radio light curve which results from these two competing effects of rapidly declining 
absorption and more slowly declining emission is shown schematically for one frequency in 
Figure 1. 
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2.1. Parameterized Radio Light Curves 

It has been shown by Weiler et al. (1986) that the radio emission from RSNe can be 
described, for simple cases, by: 



5 GHz/ V lda y/ 



where 



with Ki and K 2 corresponding, formally, to the flux density and uniform external 
absorption, respectively, at 5 GHz (6 cm wavelength) one day after the date of explosion, t . 
The term e~ T represents the attenuation of a medium that completely and uniformly covers 
the emitting source (i.e., a "uniform" external absorption) and the absorption is assumed 
due to purely thermal, ionized hydrogen with frequency dependence v~ 2A . The parameter 
S describes the time dependence of the optical depth for the uniform external absorbing 
medium. The emission from the RSN is assumed to be nonthermal synchrotron radiation 
with spectral index a and to be decreasing with time with index (5. 



2.2. Model Interpretation and Predictions 

The parameters of Equations 1 and 2 can be interpreted in terms of the Chevalier 
(1982a, b; 1984) model if the dense, external cocoon of material is established by a constant 
mass-loss rate (M), constant velocity (w) wind from a red supergiant progenitor. The time 
dependence (3 of the flux density decline is derived by assuming equipartition between the 
magnetic and relativistic particle energies, a constant ratio of those energies to the thermal 
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post-shock particle energy, a power-law density distribution for the unshocked circumstellar 
material (CSM) of Pcsm oc r~ s , and a density distribution for the unshocked ejecta of 
p e j oc r~ n . The time dependence S of the external uniform optical depth is clearly a function 
of both the shock expansion index m (R s h oc k oc t m ) and the radial distribution of the CSM 



The Chevalier model then interprets the parameters of Eq. 1 at any given frequency v 

as: 



K x oc (M/u?) (7 ~ 7+12m)/4 or K x oc {M/w) a , 



(3) 



where 



a = 



( 7 -7+12m)/4, 



(4) 



and 



P = -(7 + 5-6m)/2; 



(5) 



and of Eq. 2 as 



K 2 oc (M/w) (5 ~ 3m) or ^ 2 oc (M/w) c , 



(6) 



where 



c = 5 — 3m, 



(7) 



and 
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5 = -3m. (8) 

The term 7 is the index of the energy spectrum of the relativistic synchrotron electrons 
(which is related to the radio spectral index a by 7 = —2a +1), and m is the shock 
expansion index described above. 

Substituting these relations into Eq. 1 and including the distance to the SN to convert 
observed flux density (S v ) into spectral luminosity (L v ), we can write 

L u oc (M/w) a (t - t fexp[(M/w) c (t - t ) 5 ]. (9) 

This is a function which has a single peak when the optical depth approaches unity so that 
we can obtain a value for the mass- loss rate/presupernova wind velocity ratio (M/w) from 
setting the first derivative equal to zero (i.e., dL u /d(t — t ) = 0) which gives 

(M/w) c oc r peak = (5/5 (10) 

or 

M/w (x (P/Sy^U pcak - t )- 5 ' c . (11) 
Substituting this back into Eq. 9 yields a luminosity for frequency v at time t v pea k of: 

L v P e ak oc (M/w) a (U pcak - tofe-T" ^\ (12) 

or 
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L v peak OC (P/S) a ' C e-WQ(t v peak " * ) 



(13) 



where 



77 = — (a5 — f3c)/c. 



(14) 



Eq. 13 implies that the peak luminosity at any frequency is related to the time interval 
after explosion required for the radio emission to reach that peak. This is illustrated 
schematically in Figure 2. 

In principle, the index 77 in this simplified modelling treatment is fully determined 
by only two physical parameters, the index of the electron energy spectrum 7 (where 
7 = —2a + 1) and the index of the shock expansion m (i.e., -R s hock oc t m ). In practice, 
we do not expect to be able to describe the evolution of a SN explosion so simply and 
will employ an empirical fit to the available data. However, it is interesting to note that, 
although there is considerable scatter, the simple model prediction for the four Type Ib/c 
SNe gives, on average, relatively little variation of spectral luminosity with turn-on delay 
[L u peak oc it v P c a k — t Q ) v , i] ~ —0.2], and for the twelve Type II SNe, a much stronger 
dependence of spectral luminosity with turn-on delay (77 ~ +1.4). For the Type II SNe, it 
is quite surprising that this is the same index value as is obtained from the actual fit to 
the data in §4.2, Eq. 17. However, since the scatter in 77 values is quite large, this exact 
agreement is certainly fortuitous even though the indication of a steep slope is supported. 



2.3. Type II Peculiar SNe 



It should be noted that fitting the so-called Type II "peculiar" (Type Ilpec), or Type 
Iln, RSNe 1986J (Weiler, Panagia, & Sramek 1990) and 1988Z (Van Dyk et al. 1993b), and 
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Type lib RSN 1993J (Van Dyk et al. 1994) radio light curves requires additional absorption 
terms beyond those given by Eqs. 1 and 2. In particular, an additional term describing 
an "internal," or "clumpy external," component of the absorption, which is not included 
in the standard Chevalier model, must be used and could imply that some Type II RSNe 
must be considered separately. However, since our results are empirically determined, and 
since RSNe 1986J, 1988Z, and 1993J appear to obey the same relation as their more normal 
brethren, they are included along with other Type II SNe for the present discussion. 



2.4. SN 1987A 

One of the major peculiarities of SN 1987A is that its explosion occurred when its 
progenitor was a blue supergiant (BSG) than being a red supergiant (RSG). Because of 
the much higher wind velocity of a BSG than that of a RSG, the density of the CSM, i.e., 
the density of the presupernova stellar wind (oc M/w), around SN 1987A was considerably 
lower than around "normal" Type II SNe and, therefore, its early radio emission was 
intrinsically much weaker and evolved much more quickly than normal. In particular, the 
flux density at 6 cm reached a peak only 1 day after the explosion and fell below detection 
limits only a few days later. As a consequence, despite the efforts of Australian radio 
astronomers, the radio light curve of SN 1987A is poorly defined. 

This fact, and the large difference in circumstellar environment (much lower density, 
much higher velocity) as compared with normal Type II SNe, makes it uncertain whether a 
quantitative comparison of SN 1987A with the other Type II SNe can be made. However, 
it is apparent in Fig. 3 that SN 1987A agrees reasonably well with an extrapolation of the 
relation determined for the remaining Type II SNe which have much greater peak 6 cm 
radio luminosities and much longer times from explosion to 6 cm peak flux density. This 
may imply that the relation is indeed valid for all Type II RSNe and that the gap between 
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RSN 1996cb and RSN 1987A is due to selection effects against fast turn-on, low radio 
luminosity RSNe. Thus, since we are dealing with an empirical relation, we have included 
SN 1987A in our discussion. 

3. Data 

3.1. Peak Luminosity and Time to Peak 

In order to quantify possible relations between SN type, peak radio luminosity, and time 
from explosion to radio peak, it is necessary to establish these quantities in a systematic 
fashion. The data sets available for the various RSNe are of quite variable quality, and most 
RSNe show some deviation from smooth radio light curves. Thus, even when observations 
are available at 6 cm wavelength near the time of peak flux density, the smoother best-fit 
model values have been used in all cases to determine the two quantities of interest: peak 
6 cm flux density, Sq cm pca k (from which the peak 6 cm spectral luminosity, L 6 cm pea k, is 
calculated using the distances given in Table 1), and the time in days (t e cm pca k — to) from 
explosion, to, to 6 cm peak luminosity, t e cm pca k. These values, along with associated errors, 
are given in Table 2 for each of the 16 RSNe under consideration. 

Two caveats should be kept in mind. First, new fits have been performed for all 
available data sets and, due to revisions or the inclusion of new data, some fitting 
parameters may differ slightly from previously published values. Second, these results are 
from preliminary fits to some of the data sets and more data may somewhat alter the 
parameter values in future treatments. 



3.2. Errors 



Because of the very diverse nature of the data, missing data at critical (usually early) 
times, and variations in the assumptions for the modelling of different RSNe, it is extremely 
difficult to assign errors to the basic quantities of S 6 cm pea k, including, in many cases, large 
uncertainty in determining or assigning a date of explosion, to- Additionally, the conversion 
of Sg cm peak into L 6 cm pea k introduces the additional uncertainty of poor distance estimates 
to many of the parent galaxies of the SNe. 

However, to try to provide some indication of the relative quality of the different 
points, we have attempted to assign reasonable error estimates to all values. For the 
distances listed in Table 1, we have tried to find the best value and error estimate available 
in the literature and, when successful, have listed the appropriate reference. For galaxies 
where no estimate is available, we have taken the value from Tully (1988), accepting all his 
corrections but adjusting his distance value obtained with his assumed H = 75 km s -1 
Mpc -1 to what appears to be a more recently preferred value of H = 65 km s _1 Mpc -1 
(see, e.g., Saha, et al. 1997 and references therein; Riess, Press, & Kirshner 1996; Hamuy, 
et al. 1996). Since Tully (1988) assigns no distance error to his values, we have assumed 
that ± 15 %, to roughly cover the range of H from 55 to 75 km s _1 Mpc -1 within one 
sigma, is a reasonable estimate. 

For the other two critical parameters, S 6 cm pea k and (t 6 cm pea k - t ), we have used 
several methods to obtain rough error estimates. 

Where enough data are available and the radio light curves are well determined (e.g., 
SN 1993 J with hundreds of measurements), we were able to use the bootstrap method^ 

5 Bootstrap (Press et al. 1992) procedures use the actual data sets to generate thousands of 
synthetic data sets that have the same number of data points, but some fraction of the data is 
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to obtain direct estimates of the one sigma deviation for both quantities. However, in 
the cases where the data were very sparse (e.g., SN 1986E with only one detection at 
one frequency), we had to rely on a qualitative estimate of the error by adjusting fitting 
parameters by hand, to obtain a range of values for Sq cm pea k and (te cm pca k — to) which 
appeared consistent with the data. In extreme cases where one or both parameters were 
effectively indeterminate, that has been indicated in Table 2 by the description "ind." and, 
in Figure 3, by a line stub for that side of the range of uncertainty. 

Once error estimates were established for Sq cm pca k and the distance, determination of 
error estimates for the L 6 cm pea k followed standard error propagation procedures. 

Finally, error estimates for the "radio distances" in Table 3 were obtained for Type lb 
and Type Ic SNe from propagation of the flux density error from Table 2 and the deviation 
from the average luminosity (Equation 15 for Type lb and Equation 16 for Type Ic) for 
each SN type. For Type II SNe errors were estimated from the flux density error from Table 
2 and the standard deviation of the luminosities from the straight line shown in Figure 3. 
While it is felt that these error estimates are reasonable for the quality of the data available, 
they cannot be rigorously justified. 



4. Discussion 



replaced by duplicated original points. The fitting parameters are then estimated for these 
synthetic data sets using the same algorithms that are used to determine the parameters 
from the actual data. The ensemble of parameter fits is then used to estimate errors for the 
parameters by examining the number distribution of the parameter in question. The errors 
we present correspond to the 15.85% and 84.15% points in the distribution of the occurrence 
of parameters from the synthetic data. 
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4.1. Type Ib/c Peak Radio Luminosities and Distances 

The model discussion presented in Section 2 suggests that Type Ib/c SNe may have 
relatively little variation in their peak spectral luminosities, even if there are variations in 
the measured time from explosion to peak. Unfortunately, there are only four Type Ib/c 
RSNe which have been studied at radio wavelengths, and even these four examples are 
diverse, with two being Type lb and two being Type Ic. Additionally, SN 1984L has very 
limited data available, and SN 1990B was quite distant and therefore studied with poor 
signal-to- noise. Within these limitations, however, examination of Table 2 shows that Type 
lb RSNe, based on data for only two objects, may be approximate standard candles with 
peak 6 cm luminosities ranging only from 14.1 - 25.7 x 10 26 erg s" 1 Hz" 1 , and, similarly, 
Type Ic RSNe, again based on data for only two objects, may be approximate standard 
candles with peak 6 cm luminosities ranging only from 5.6 - 13.7 x 10 26 erg s' 1 Hz -1 . 
These are consistent to within the errors for Type lb SNe having an average peak 6 cm 
luminosity of 

U cm P eak « 19.9 x 10 26 erg s" 1 Hz" 1 (15) 
and Type Ic SNe having an average peak 6 cm luminosity of 

L 6 cm peak ~ 6.5 x 10 26 erg s" 1 Hz" 1 . (16) 

If we assume that Type lb and Type Ic RSNe are, in fact, standard candles in their 
peak 6 cm radio emission and that, on average, the assumed distances listed for the objects 
in Table 1 are correct, we can determine a "radio distance" to each RSN, such that its 
measured peak 6 cm flux density would yield the average peak 6 cm spectral luminosity for 
its particular type. These distances are listed in Table 3 along with an error. The "assumed 
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distances" and errors from Table 1 are repeated in Table 3 for the reader's convenience. 

While it should be kept in mind that the two sets of distance estimates are not 
completely independent, with the assumed "independent" distances de facto setting the 
scale for the radio distances, it is apparent that the data are consistent, to within the errors, 
with an assumption of constant 6 cm peak luminosity. Thus, the distance to any Type 
lb or Type Ic RSN can be estimated simply by measuring its 6 cm peak flux density and 
comparing this observed value with the luminosity predicted by Eqs. 15 and 16. Clearly, 
more objects of both subtypes must be identified and studied in the radio to better test and 
define these two equations. 



4.2. Type II Peak Radio Luminosities and Distances 

There are 12 examples of Type II RSNe which have a sufficient number of radio 
observations available for detailed study. While Type II RSNe are far more heterogeneous 
in their radio (as well as in their optical) properties than Type Ib/c RSNe, inspection of 
Table 2 and Figure 3 shows that the longer it takes a Type II RSN to reach 6 cm peak flux 
density, the higher its 6 cm luminosity at that peak. This is also supported by the model 
discussion in Section 2, which implies a relatively steep power-law relation between the two 
quantities. However, as with the Type Ib/c SNe, we are considering here an empirical test 
and therefore will determine a best-fit power-law to the available turn-on time vs. peak 6 
cm luminosity data. 

Performing an unweighted least squares fit to the data for the 12 available Type II 
RSNe, we obtain the relation (formal fitting errors given) of 

L 6 cm peak ^ 5 - 5 (^3 4) x 1Q23 ( t6 cm P- k " *o) 1 - 4±0 - 2 erg s" 1 Hz" 1 , (17) 
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with time in days. This relation is shown as the dashed line in Figure 3. We note that the 
empirically derived slope of the L 6 cm pcak vs. (t e cm pcak — t ) relationship agrees well with 
the value inferred from model estimates (see Eq. 14). This result suggests the validity of 
the Chevalier (1982a,b) model, at least for interpreting the origin of the radio emission. 

If we assume that Type II RSNe obey this relation and that, on average, the 
independent distances listed for the objects in Table 1 are correct, we can determine a 
"radio" distance to each RSN, such that its measured peak 6 cm flux density would yield 
the "best-fit" 6 cm spectral luminosity predicted by Eq. 17. These "radio distances" are 
listed in Table 3, Column 4, along with errors, while the "assumed distances" and errors 
from Table 1 are repeated in Table 3 for the reader's convenience. 

Again, while it should be kept in mind that the two sets of distance estimates are not 
completely independent, with the "assumed distances" de facto setting the scale for the 
radio distances, it is apparent that the data are consistent, to within the errors, with the 
results of Eq. 17. Thus, the distance to any Type II RSN can be estimated simply by 
measuring its time from explosion to 6 cm peak flux density and comparing the observed 6 
cm peak flux density with the predicted luminosity. 

5. Deviations 

Examination of the distances determined from the relations of Eqs. 15, 16, and 17 for 
Type lb, Type Ic and Type II, respectively, in Table 3 shows that they are generally in 
good agreement to within the estimated errors. However, several cases show relatively large 
disagreement inconsistent with the error estimates. 
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5.1. SN 1985L, SN 1986E, and SN 1996cb 

Perhaps too easily, large deviations for these three SNe can be explained, if not 
dismissed. SN 1985L was faint (S 6 cm pca k ~ 0.65 mJy, estimated) and had only two clear 
detections of radio emission (Van Dyk et al., 1998). While reasonable estimates of RSN 
parameters can be obtained for SN 1985L through the use of a number of upper limits, the 
quality of the L 6 cm pca k and (te cm pca k — t ) must remain suspect. 

SN 1986E (Montes et al. 1997) has the same difficulties as SN 1985L, but even more 
extreme. Only one clear radio detection at one frequency exists, and the rest of the 
parameter fitting has to rely on the use of upper limits. It was also quite faint, with a peak 
6 cm flux density of only Sg C m peak ~ 

0.31 mJy. 

SN 1996cb (Van Dyk et al., in preparation), while relatively bright with a 6 cm peak 
flux density of S6 cm peak ~ 1.8 mJy, exhibited the most rapid turn-on of any Type II RSNe, 
except for SN 1987A. (In fact, many of its radio properties resemble those of the Type Ib/c 
RSNe.) This, unfortunately, meant that the turn-on at wavelengths shorter than 20 cm was 
completely missed, and estimates of S 6 cm pea k and (t 6 cm pca k — t ) are rather uncertain. 

5.2. SN 1980K in NGC 6946 

RSN 1980K (Weiler et al. 1986, 1992) has a very large deviation from Eq. 17 (the 
dashed line in Figure 3). Also, it is possibly the most interesting deviation in that it is not 
easily explained away. SN 1980K was a reasonably strong, well-studied, "normal" Type II 
RSN, where there can be little doubt as to its characteristics and model parameters. Thus, 
it may well serve as a test case. 

If the relation of Eq. 17 for Type II SNe is valid, and SN 1980K is indeed a normal 
Type II RSN, then the distance estimate of 6.3 Mpc from Tully (1988; converted to H = 
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65 km s 1 Mpc 1 ) for the parent galaxy, NGC 6946, must be too short. A significantly 
larger distance of ^ 10 Mpc is implied by the radio results. 

NGC 6946 is a large (11.5 x 9.8 arcmin), nearby (^heliocentric = 46 km s" 1 ; Tully 1988) 
spiral (Scd) galaxy, which has been a prolific producer of SNe in the last century (SNe 
1917A, 1939C, 1948B, 1968D, 1969P, and 1980K). However, it is at a low galactic latitude 
(b = 11.7°) and located along the direction of rotation of the Milky Way (I = 95.7°), 
requiring large corrections to obtain an intrinsic radial velocity with respect to the galactic 
rest frame. Because it is so close, peculiar motions within the Local Group of galaxies may 
also disturb the Hubble flow and yield an inaccurate distance estimate from corrected radial 
velocity measurements. 

Recent work on the Expanding Photosphere Method (EPM) by Schmidt, Kirshner, & 
Eastman (1992) yields a distance to SN 1980K (and thus, NGC 6946) of 8.1 ± 1.5 Mpc from 
IR observations and 7.2 (+0.7, —1.0) Mpc from optical observations, both somewhat closer 
than our radio estimate of 12.3 ±1.9 Mpc. Furthermore, a revision of the EPM results by 
Schmidt et al. (1994) yields a distance of 5.7 ± 0.7 Mpc, in much closer agreement with the 
Tully (1988) distance than with our radio distance. 

We have circumstantial evidence that NGC 6946 may be more distant than commonly 
assumed. The fact that NGC 6946 is the galaxy with the largest number of SN events 
ever recorded suggests a high blue luminosity, which should be at least as high as that of 
another prolific SN producer, NGC 4321 (M100), in which 4 SNe have been discovered 
(with only one, SN 1979C, being a bona fide Type II SN). If we assume that NGC 6946 
has the same blue luminosity, M B) as NGC 4321 and establish the difference in their total 
apparent magnitudes as m B (NGC 4321) - m B (NGC 6946) = 0.47 (Tully 1988), accepting 
that the difference of their reddening corrections is ^1 magnitude (as indicated by optical 
and UV observations of SNe 1979C and 1980K; Panagia 1982), we find that the distance 
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to NGC 6946 is no less than 0.51 times the distance to NGC 4321 (d = 17.1 ± 1.3 Mpc; 
Freedman et al. 1994a). Thus, this implies a distance to NGC 6946 of d > 8.7 Mpc. 

In any case, it appears that the distance to NGC 6946 is sufficiently uncertain that it 
may be in error by the factor of ~ 2 which the radio data for SN 1980K imply. Clearly, a 
measurement of the distance to NGC 6946 with more accurate techniques (e.g., Cepheids) 
is called for and could provide an important test of the L 6 cm pca k vs. (te cm pC ak — to) 
hypothesis. 

6. Conclusions 

We have presented evidence that the radio emission from SNe may have quantifiable 
properties which allow for distance determinations. Type lb RSNe, based on a statistically 
very small sample of only two objects, may be approximate radio "standard candles," with 
with a 6 cm peak luminosity of ~ 19.9 x 10 26 erg s -1 Hz^ 1 ; Type Ic RSNe, also based on a 
very small sample of only two objects, may be approximate radio "standard candles" with 
a 6 cm peak luminosity of ~ 6.5 x 10 26 erg s -1 Hz -1 ; and Type II RSNe, based on a sample 
of twelve objects, appear to obey a relation L 6 cm pca k — 5.5 x 10 23 (te cm peak — ^o) 1 ' 4 erg 
s _1 Hz -1 (with time in days). Thus, measurement of the radio turn-on time (t 6 cm pea k — 1 ) 
and peak flux density S e cm pca k can yield a luminosity estimate and therefore a distance. 

The reality of the "standard radio candle" hypothesis for Type lb and Type Ic RSNe 
may be tested simply through the study of more objects. For Type II RSNe, the large 
deviation of the well-observed SN 1980K in NGC 6946 from the distance predictions of this 
hypothesis could serve as a useful test. Whereas the presently accepted distance is ~6 Mpc 
to NGC 6946 (Tully 1988), the radio results predict a significantly larger distance of ^10 
Mpc. A new, accurate determination of the distance to NGC 6946 is certainly called for 
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and, in general, improved distances to all host galaxies of RSNe would be valuable. 

Although there are still relatively few objects to which these techniques can be applied, 
RSNe may eventually provide a powerful and independent technique for investigating the 
long-standing problem of distance estimates in astronomy. With such intrinsically bright 
Type II RSNe as SN 1988Z and SN 1986J, the technique can be applied to distances of 
at least 100 Mpc with current VLA technology. With future sensitivity improvements 
and planned, new, more sensitive radio telescopes, the technique can be extended to large 
distances even for less luminous RSNe. With a sensitivity of 1 /iJy, a SN of the same class 
as SNe 1988Z and 1986J could be detected up to a redshift of ~ 1, while a more normal 
Type II SN, such as SN 1980K, could be studied accurately up to a redshift of ^ 0.1. 



KWW & MJM wish to thank the Office of Naval Research (ONR) for the 6.1 funding 
supporting this research. Additional information and data on radio supernovae can be 
found on [http: / / rsd- www. nrl. navy.mil / 72 14/ weiler /| and linked pages. 
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Table 1. Radio Supernova List and Parent Galaxy Properties 

RSN Optical Parent Hubble Systemic Assumed Dist. 

Name Type Galaxy Type Velocity a Distance 13 Ref. c 

(km s- 1 ) (Mpc) 



Type Ib/c 



1983N 


lb 


NGC 


5236 


(M83) 


ScI-II 


337 


5.4 ±0.8 


1984L 


lb 


NGC 


991 




Sell 


1541 


21.7 ±3.3 


1990B 


Ic 


NGC 


4568 




Sbc 


2168 


19.4 ±2.9 


19941 


Ic 


NGC 


5194 


(M51) 


Sbc 


573 


8.9 ±1.3 


rype II 
















1970G 


IIL 


NGC 


5457 


(M101) 


Scd 


379 


7.4 ±0.6 1 


1978K 


II 


NGC 


1313 




Sd 


254 


4.5 ±0.7 


1979C 


IIL 


NGC 


4321 


(M100) 


Sbc 


1522 


17.1 ±1.8 2 


1980K 


IIL 


NGC 


6946 




Scd 


338 


6.3 ±1.0 


1981K 


II 


NGC 


4258 


(M106) 


Sbc 


521 


5.9 ±1.1 3 


1985L 


IIL 


NGC 


5033 




Sc 


931 


21.6 ±3.2 


1986E 


IIL 


NGC 


4302 




Sc 


1044 


19.4 ±2.9 


1986J 


Iln 


NGC 


891 




Sb 


712 


11.1 ± 1.7 


1987A 


Ilpec 


LMC 






Irr 


3.9 


0.051 ±0.002 4 



Table 1 — Continued 



RSN Optical 
Name Type 



1988Z Iln 
1993J lib 
1996cb lib 



Parent 
Galaxy 



MCG +03-28-022 
NGC 3031 (M81) 
NGC 3510 



Hubble Systemic 
Type Velocity 21 
(km s" 1 ) 

S 6670 
Sab 96 
SBm 770 



Assumed Dist. 
Distance 13 Ref. c 
(Mpc) 



102.6 ±15.4 5 
3.6 ±0.3 6 
9.1 ± 1.4 



a From Tully (1988); heliocentric velocity of the parent galaxy is corrected for Galactic 
rotation of 300 km s" 1 towards I = 90°, b = 0°. 

b Unless otherwise referenced, the distance is from Tully (1988), and is based on 

systemic velocity, corrected by a model which assumes the Milky Way is retarded by 300 

km s -1 from universal expansion by the mass of the Virgo Cluster, and on an assumed 

Hubble constant of H = 65 km s" 1 Mpc" 1 [Tully (1988) uses H = 75 km s" 1 Mpc" 1 ]. 

Unless otherwise referenced, errors are conservatively taken to be ±15%, a value which 

is large enough to encompass both H = 55 and H = 75 km s _1 Mpc -1 within 1 a. 
c (l) Kelson et al. 1996; (2) Freedman et al. (1994a); (3) average of VLBI H 2 maser 

dynamical distances from Greenhill et al. (1995) [5.4 ± 1.3 Mpc] and from Miyoshi et 

al. (1995) [6.4 ± 0.9 Mpc]; (4) Panagia et al. 1991; (5) systemic velocity from Stathakis 

& Sadler (1991), and distance derived from this velocity and an assumed H Q = 65 km 

s- 1 Mpc" 1 ; (6) Freedman et al. (1994b). 
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Table 2. Radio Supernova Properties 



RSN Opt. Peak 6 cm Time from Explosion Peak 6 cm RSN 

Name Type Flux Density a to 6 cm Peak a Luminosity a ' b Refs. 

('Sfj cm peak) (^6 cm peak ^o) (-^6 cm peak) 

(mJy) (days) (10 26 erg s" 1 Hz" 1 ) 



Type Ib/c 

1983N lb 

1984L lb 

1990B Ic 

19941 Ic 



Type II 

1970G 

1978K 

1979C 

1980K 

1981K 

1985L 

1986E 

1986J 

1987A 



IIL 

II 

IIL 

IIL 

II 

IIL 

IIL 

Iln 

Ilpec 



40.1 (3?) 
4.6 ( + -'f„) 

1-3 (3!) 
14.6 (3J) 



2i.6 (3:1) 

229.0 («I;J) 
7-3 

2.5 (31) 
5.1 («) 

o.7 (3:1) 
0.3 (31) 

136.0 (!«) 
91.4 ( + _2Zf) 



H.6 (3J) 

n.o (Sf) 
37.5 ft-) 

35.8 (3;») 



307.0 («J 



685.0 («$ 



605.0 (*J 



140.0 



-29.3 
-27.7 



33-8 (? n l) 
309.0 ( + _ 1 f ° 



224.0 
1150.0 



-32.0 

+85.0 
-21.0 



14-1 (33 

25-7 (3£) 
5.6 (33 
13.7 (!») 

14.0 (3J) 

55.3 (31;) 

25.3 
1.2 



+5.4 
-5.4 



2.1 



-0.4 
-0.4 

/+ind 



findA 
-1.2 J 

3.6 

199.0 

0-003 ft;™) 



1,2 
1,3 

4 

5 

1,6,7,8 

9,10 
1,11,12 

1,13 
1,14,15 
16 
17 
18 
19 
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Table 2 — Continued 



RSN Opt. Peak 6 cm Time from Explosion Peak 6 cm RSN 

Name Type Flux Density a to 6 cm Peak a Luminosity a ' b Refs. 

cm peak ) (^6 cm peak ~~ ^o) (L 6 cm peak ) 

(mJy) (days) (10 26 erg s" 1 Hz" 1 ) 



1988Z Iln 1.9 (tJJ) 1420.0 (i™-° ) 237.0 20 

1993J lib 95.2 180.0 15.0 21 

1996cb lib 1.85 (+J; 2 ) 19.4 (+^) 1.8 (+^) 22,23 



a See §3.2 for discussion of errors. 

b Distances are from Table 1 and peak flux densities are from Column 3. 

c (l) Weiler et al. (1986); (2) Sramek, Panagia, & Weiler (1984); (3) Panagia, Sramek, & Weiler 
(1986); (4) Van Dyk et al. (1993a); (5) Rupen et al. (1998); (6) Allen et al. (1976); (7) Marscher & 
Brown (1978); (8) Cowan, Goss, k Sramek (1991); (9) Ryder et al. (1993); (10) Montes, Weiler, & 
Panagia (1997); (11) Weiler et al. (1981); (12) Weiler et al. (1991); (13) Weiler et al. (1992); (14) 
van der Hulst et al. (1983); (15) Van Dyk et al. (1992); (16) Van Dyk et al. 1998; (17) Montes 
et al. (1997); (18) Weiler, Panagia, & Sramek (1990); (19) Turtle et al. (1987); (20) Van Dyk et 
al. (1993b); (21) Van Dyk et al. (1994); (22) Van Dyk et al. (1996); (23) K. W. Weiler, private 
communication. 
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Table 3. Distances 



RSN Optical Assumed Radio 

Name Type Distance 51 Distance 13 

(Mpc) (Mpc) 



Type Ib/c 



1983N 


lb 


5.4 


± 


0.8 


1984L 


lb 


21.7 


± 


3.3 


1990B 


Ic 


19.4 


± 


2.9 


19941 


Ic 


8.9 


± 


1.3 



6.5 

i9.i 



+3.2 
-3.3 



20.7 

s i (US) 

Type II 

1970G IIL 7.4 ± 0.6 7.4 

1978K II 4.5 ± 0.7 3.9 (+™) 

1979C IIL 17.1 ± 1.8 20.1 

1980K IIL 6.3 ± 1.0 12.7 

1981K II 5.9 ± 1.1 3.3 (t^J) 



1985L IIL 21.6 ± 3.2 42.5 

1986E IIL 19.4 ± 2.9 49.0 



+7.2 
-9.8 

+8.3 
-8.3 



1986J Iln 11.1 ± 1.7 7.2 

1987A Ilpec 0.051 ± 0.002 0.071 
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Table 3 — Continued 



RSN Optical Assumed Radio 

Name Type Distance 21 Distance 13 

(Mpc) (Mpc) 



1988Z Iln 102.6 ± 15.4 70.7 

1993J lib 3.6 ± 0.3 2.4 (+^) 

1996cb lib 9.1 ± 1.4 3.8 



a From Table 1. 

b Distance which the RSN would have if: 1) for Type Ib/c RSNe, the peak observed 6 
cm flux density originates from a "standard candle" radio source with peak luminosity 
given, respectively, by Eqs. 15 and 16; or 2) for Type II RSNe, the peak observed 6 cm 
flux density originates from a source with peak luminosity given by Eq. 17, with time in 
days. 
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Fig. 1. — Schematic representation of the competing effects of rapidly declining absorption 
(long-dashed curve) and more slowly declining emission (short-dashed curve), which yield the 
characteristic rapid turn-on, slower turn-off "light curve" (solid curve) observed for RSNe at 
any given frequency. 



-30- 



Fig. 2. — Schematic representation of the radio light curves at a given frequency for SNe 
with differing mass- loss rate/presupernova stellar wind velocity (M/w) ratios. Note that the 
longer the time delay required to reach peak luminosity the higher is that luminosity peak. 
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Fig. 3. — Peak 6 cm luminosity, L 6 cm peak of RSNe vs. time, in days, from explosion to peak 
6 cm flux density (te cm peak — to). Type II SNe are plotted as filled triangles. The dashed 
line is the unweighted, best fit to the 12 available Type II RSNe. The Type II SNe show a 
large range in times to 6 cm peak (t 6 cm pcak — t ) and in peak 6 cm luminosity L 6 cm pcak , 
but appear to obey the relation given by Eq. 17. Error bars are based on best estimates 
(see §3.2). Where no error or only a stub of a line is shown, the error in that direction is 
indeterminate. 
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